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Abstract 

We discuss full-sky convolution of the instrumental beam with the CMB sky signal in CMB 
single-dish and interferometry experiments, using the method 13] that the measured temperature 
and polarization anisotropics are defined globally on the group manifold of the three-dimensional 
rotation by means of Wigner D-functions. We re-derive the anisotropy and polarization correlation 
functions incorporated with asymmetric window functions, which are then explicitly calculated for 
a single-dish elliptical Gaussian beam and an interferometric Gaussian beam. 
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I. INTRODUCTION 



Since the detection of the large-angle temperature anisotropy of the cosmic microwave 
background (CMB) by the Cosmic Background Explorer (COBE) satellite 0, many CMB 
measurements have reported detections or upper limits of the CMB temperature anisotropy 
power spectrum Cri over a wide range of Z . Recently, the Wilkinson Microwave Anisotropy 
Probe (WMAP) [3| has measured the Cti in an unprecedented accuracy for / < 800 and 
thus made a precise estimation of a number of cosmological parameters. 

Polarization of the CMB contains a wealth of information about the early universe as 
well. It can cross-check the measured Cti and improve the accuracy in determining the 
cosmological parameters. A CMB polarization field can be decomposed into an electric-type 
E mode and a magnetic-type B mode Recently, the DASI instrument, a ground-based 
interferometric array with degree-scale resolution, has detected the CMB i^polarization and 
TE cross-correlation, while setting an upper limit on the B polarization [5|. Furthermore, 
the WMAP has measured the TE power spectrum 0, which is consistent with theoretical 
predictions based on the measured CMB anisotropy and indicates a significant large-scale 
E polarization. Measuring the CMB polarization has become one of the main goals of 
CMB experiments . Future ground-based and balloon-borne experiments and the Planck 
satellite will unveil detailed features of the CMB anisotropies, thus allowing one to determine 
to a high precision the cosmological parameters. 

As high-resolution data made by high-sensitivity CMB experiments comes along, it will 
be necessary to consider the effect of beam asymmetry in data analysis in order to make 
unbiased estimation of the power spectra and ultimately of the cosmological parameters. 
This effect, mainly caused by the off- axis position of the detector in the focal plane of the 
telescope, is small and so far has been largely neglected in CMB data analysis. Recently, the 
effect of beam asymmetry has been investigated. Numerical studies III have shown 

that the typical difference due to beam asymmetry is of few fiK in Planck configuration and 
increases with the beamwidth and ellipticity 0] and that an azimuthally symmetrized beam 
does not bias the power spectrum estimates for most practical situations 0. In Ref. 
a fast convolution algorithm based on the decomposition of the asymmetric beam as a sum 
of circular Gaussian functions is developed. On the other hand, analytic methods have 
been formulated for treating arbitrary asymmetric beams in full-sky CMB temperature and 
polarization anisotropy experiments by using spin-weighted spherical harmonics 2| and 



Wigner D-functions 13]. In Ref. |13[, the authors have pointed out that for asymmetric beam 
functions the result of convolutions is on the group manifold of the three-dimensional rotation 
and provided an algorithm for fast computation of the convolution on this group manifold. 
This method has been applied to study the effect of an elliptical Gaussian beam in the 
estimation of temperature and polarization correlation functions using semi-analytic or full 
numerical integration [l^ and a perturbative series in powers of the ellipticity parameter [15| . 

In this paper, we will study the effect of beam asymmetry in full-sky CMB single-dish 
as well as interferometry experiments. We will derive the full-sky convolution for a com- 
plex asymmetric beam and define the measured temperature and polarization anisotropies 
globally on the group manifold of the three-dimensional rotation by using Wigner D- 
functions . As such, we can extend the standard anisotropy and polarization correlation 
functions of domain over the celestial sphere to those over the rotation group manifold. The 
result is then applied to calculate the full-sky correlation functions for two typical cases: 
an elliptical Gaussian beam and an interferometric Gaussian beam. The former is a good 
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approximation to the actual shape of the window in most of CMB single-dish experiments 
and the latter is maximally asymmetric due to the finite length of the baseline of the inter- 
ferometer. We will be able to obtain closed forms for the covariance matrices in both cases 
which are directly applicable in the CMB likelihood data analysis. 

The paper is organized as follows. In Sec. IHl a brief account of the building block of CMB 
anisotropics is given. We present the full-sky convolution of a complex asymmetric beam 
with the CMB sky in Sec. 11111 The full-sky CMB anisotropy and polarization correlation 
functions are re-derived in Sec. IIVI In Sec. El we apply the result to study the effect of 
beam asymmetry in two cases: a single-dish elliptical Gaussian beam and an interferometric 
Gaussian beam. The flat-sky limits of these two cases are discussed in Sec. IVII Sec. I VI II is 
our conclusions. 



II. CMB TEMPERATURE AND POLARIZATION ANISOTROPIES 

Polarized emission is conventionally described in terms of the four Stokes parameters 
(/, Q, U, V), where I is the intensity, Q and U represent the linear polarization, and V de- 
scribes the circular polarization. Since circular polarization cannot be generated by Thomson 
scattering alone, the parameter V decouples from the other components and will not be con- 
sidered. Let us define T be the temperature fluctuation about the mean; then, the CMB 
anisotropics are completely described by {T,Q,U), where each parameter is a function of 
the pointing direction e{6, (p) on the celestial sphere. 

Considering the CMB as Gaussian random fields, we can expand the Stokes parameters 

as 



T{e) = '^ao^imYiraie), 

Im 

{Q-iU){e) = 2>^m(e), 

Im 

{Q + tU){e) = Y.('-2M-2Yim{e), (1) 

Im 

where ao,im and a+o . /rp, are Gaussian random variables, and ±2^Zm are spin-2 spherical har- 

2/ + l(/ + m)!(/-m) 



monies given by [16[ [25[ 



sYlm{0, ' 





Att {l + s)\ {l-s)\ 

-l)'-'^-"cot'''+'-™ ( - ) , (2) 



r + s — m 



where 



max(0, m — s) < r < min(Z — s,l + m). (3) 
Note that sYim is related to Wigner D-functions Dl^,^ [l^ by 



DL^mii'. 0, 0) = ''^4,^(^)e— ^ where ,¥^9, 0) = (-1)^+'" J-^rf^^(^). (4) 
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It is straightforward to show the conjugation and symmetry relations 



(5) 
(6) 



Isotropy in the mean guarantees the ensemble averages: 



(^2,l'm'(^2,lm/ — {C El + C Bl)^Vl^m' mi 

^0,l'm,'^2,lm} ~ ^C^ciSinSm'm, (7) 

where Cti, Cei, Cbi, and Cci are respectively the anisotropy, E polarization, B polarization, 
and TE cross correlation angular power spectra. 

Consider two pointings e and e' on the celestial sphere. By making use of Eq. ((Tj) and 
the generalized addition theorem 2| 



'2/ + 1 
47r 



(8) 



we find the two-point correlation functions [2 

2/ -I- 1 

{T*{e')T{e)) =j:^CTiPi{cosf3), 



(T*(e') [Qie)+^Uie)]) = -Y: 



21 + 1 



f 47r \ (/ + 2 



(/-2)! 



-2ia 



[Q(e')+^Uie')nQie)+^Uie)]) = Y: 
[Q{e')~^U{e')r[Q{e)+^U{e)]) = Y: 



121 + 1 

l2l + 1 
An 



{CEi + CBi)2Yi-2{P,0)e- 



-2j(a— 7) 



{Cei -Cbi) 2>^2(/9,0)e 



-2i(a+7) 



(9) 



where a, and 7 are the angles defined in Fig. d Therefore, the statistics of the CMB 
anisotropy and polarization is fully described by the four independent power spectra or their 
corresponding correlation functions. The details about the evaluation of the power spectra 
can be found in Ref. 



III. FULL-SKY COMPLEX ASYMMETRIC BEAM CONVOLUTION 

In realistic CMB observations, as a result of the finite beam size of the antenna and the 
beam switching mechanism, a measurement is actually a convolution of the antenna response 
with the CMB Stokes parameters, 

Xf'^P(e,M)= J de'R{e';e,u)Xs{e'), (10) 
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FIG. 1: Spherical coordinates showing two unit vectors e'{9',(p') and e{6,(j)) with separation angle 
p. The angles between the great arc connecting the two points and the longitudes are 7 and a. u 
is an arbitrary reference axis on the asymmetric beam pattern or an interferometric baseline vector 
tangential to the sphere with the orientation angle ip. 



where Xq = T, X±2 = Qt^U, and R{e'; e, u) denotes the response function with the pointing 
direction e and the orientation direction u {u = uu, see Fig. Q). The orientation angle is 
denoted by ip and therefore we have O<0<7r, O<0< 2tt, and < ip < 2vr. Note that 
R{e'; e, ii) may be an arbitrary complex function. Expanding 



sYimie'A') 



I A-rr 

E sYim'{(3, a) e^^^ -m'Yim{e, 0) (11) 



in Eq. (Q), we obtain the full-sky convolution for an arbitrary beam function as 2 



Imm' 



In typical CMB measurements, the field of view is small such that we can make a local fiat- 
sky approximation = 7. After changing the integration variable a ^ a — %p m Eq. ()12|l 
and using Eq. (jHl), the orientation angle ip of the response function is absorbed in the new 
variable and we thus obtain 



XrP(^,0,^) = e^^"^ E asMDtm'{<P.O,i,)J'^^ Am', (13) 

Imm' 



skm' = J I sin P dj3da R{(3, a) sYim' {P, a) 



This equation was first derived in Ref. [ij by using an alternative way. From the conjugation 
relation (0), we have sbim = (—1)^"'"™' ^sh-m- We define a spin-s window function 

,Wi = Y.\Mm?. (14) 
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For a simple single-dish experiment with an axisymmetric Gaussian response function given 

by 

where cxb -C 1 is the Gaussian beamwidth, it can be shown from Eq. (O that |il Q 

.C= sW^h.,^,,, (16) 

XrP(^,0) = T^^sMsW^'^ sYlm{9,(^), (17) 
Im 



where sWi^ = exp {- [/(/ + 1) - s^] a^} . (18) 



For s = 0, ±2 and cr^ 0, the Gaussian window function gWl^ 1 and Eq. (fT7|) reduces 
to Eq. dH). 

IV. CORRELATION FUNCTIONS ON ROTATION GROUP MANIFOLD 

The result in Eq. p3p shows that CMB measurements with asymmetric beam, although 
the sky is approximated as flat locally, can be interconnected globally on the celestial sphere 
extended with a circle space of the orientation angle of the asymmetric beam. More precisely, 
the observed anisotropics are well-defined on the group manifold of the three-dimensional 
rotation Let us denote a point on the manifold by r = (e, u) = (6*, 0, ■?/'), then by 

making use of Eqs. ((Tj), ©, and (fT5|) . we obtain the two-point correlation functions (see also 
Refs. ill. Hi ) 



(X7P*(f')XrP(f)) = e-'^'^V^^ Yl J^^C7,,,/-l)™'+™„,F,_„(/5,0) X 

Im'm ' 

e-^"^(-^)e'"''^^-^"> s'bL, skm, (19) 

where Coo,« = Ct,i, Co-2,1 = —Cc,i, and C^2-2,i = Ce,i ± Cb,i- For a Gaussian beam given 
in Eq. (fT^ . substituting the result from Eq. for both s'bim' and gbim in Eq. (fT^. we 
obtain 



(Xr^ *(r-")XrP(f)) = V ^C'.,,,, ^s'MP, 0)e-"e-'^ s'W^'^ ^W^^ (20) 

where sWf' is the Gaussian window function given in Eq. (jl8|) . This reduces to the correlation 
functions in Eq. Q as the Gaussian beamwidth o";, ^ 0. 

To write Eq. (fT!I|) in a compact form, we define Xf^^^{r) = e'"^^^ Xf^^^ (r) , a = a — ip, and 
7 = 7 ~ X where ip' = it + x- After using Eq. Q, we obtain 

(xr^ *(r-")Xr'^(r)) = E -^C,s,i{-ir'DLm'{c^,f3, -7) s'&L' .fe/r.. (21) 

Im'm 

Note that a and 7 are respectively the angles between the orientation directions at the 
pointings e and e' and the great arc connecting the two pointing directions (see Fig. Q]). 
Furthermore, the root-mean-square anisotropy and polarization fluctuations are given by 
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the correlation functions ()21|) at zero lag. By taking r'^r and making use of the limiting 
property (Eq. (4.16.2) of Ref. Il7[) 

DL^,(a,0,7) = 5we-^"^("+^), (22) 

we have a — 7 = — tt and 

(i^r^(r)r) = E^^^^M.^. (23) 

where is given by Eq. (fT^ . 

Eq. (j?T|) is useful for constructing the covariance matrices in the likelihood analysis of 
CMB data made by an asymmetric beam. In the next section, we will show that in most 
practical cases the summation in Eq. pijl can be largely reduced to a tractable one. For 
examples, in the case of a slightly elliptical Gaussian beam, the summation over m for a 
fixed / converges very fast as | sbim\l\ sbio\ falls off rapidly with increasing \m\. In the case of 
an interferometric beam with a long baseline | sbim\ is approximately dependent on I only. 



V. ASYMMETRIC WINDOW FUNCTIONS: TWO CASES 

In deriving the main result (|13|). we have approximated the sky as locally flat. This is 
a good approximation for CMB experiments with high spatial resolution or small fields of 
view. We can further evaluate the integral for sbim in Eq- fjl^p by use of the approximation 
(Eq. (4.18.1.2) of Ref. [13), 

ci^„(e) ^ (-l)^+'"J.+™(/e) for /»1, (24) 

where J„(x) are Bessel functions and we have used the property J-n{x) = (— l)"J„(x). 
Below we discuss two typical cases in CMB experiments. 



A. Single-dish Elliptical Gaussian Beam 



For most of simple CMB single-dish experiments, the response function can be approxi- 
mated by an elliptical Gaussian beam. 



R^'^iP, a) 



1 



X 



2'7iax(Jy 



exp 



2al 



2al 



(25) 



where x = (3 cos a, y = (3 sin a, and and ay are the beamwidths in the major and 
minor axes respectively. Here we assume cr^ <^ 1 and -C 1. Substituting the beam 
in Eq. (jl3p . making the approximation (j24|) in the integral, and using Eqs. (3.915.2) and 
(6.651.6) of Ref. jl8| to perform the a-integration and /^-integration respectively, we obtain 



1 

2 

for 



1 + (-1)^+^ 
s + m > 



exp 



1 and 1 > (a^ 



(26) 



where In{x) are modified Bessel functions. For s = 0, the result coincides with that for 
an unpolarized elliptical Gaussian beam found in Ref. jl9|. In the limit of ax 



ay, smce 
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/(s+m)/2(0) = (5_m,s, we obtaiii that sbj^ — exp{—Pay/2)S-m,s- This is expected for an 
Gaussian beam as discussed in Eq. ()16|). From Eq. ()26|) . the eUiptical Gaussian window 
function is given by 



EG 



EG 



exp 



exp 



[1/2] 

ll{z) + 2Y.ll{z) 

n=l 

[l/2\-l 



ll{z)+2 ^n(^)+^[V2](^)+/[Wl 



n=l 



(27) 

(28) 



where z = ^^(cr^ — (Ty)/4 and [1/2] denotes the integer part oi 1/2. Making use of the relation 
(Eq. (8.537.2) of Ref. ^) 



J2 Jk{z)Jn-k{z) = Jni'^z) and Ik{z) = i ''Jkiiz), 

k=—oo 



(29) 



we find that for 2; < 1, 



- exp [-l\al + al)/2\ I, {l\al - aJ)/2) . 



(30) 



B. Interferometric Gaussian Beam 

Recent advancement in low-noise, broadband, GHz amplifiers has made interferometry a 
particularly attractive technique for detecting CMB anisotropies. An interferometric array 
is intrinsically a high- resolution polarimetric instrument, well suited to observing small-scale 
polarized intensity fluctuations while being flexible in the coverage of a wide range of angular 
scales, with resolution and sensitivity determined by the aperture of each element of the array 
and the baselines formed by the array elements. In addition, many systematic problems 
inherent in single-dish experiments, such as ground and near field atmospheric pickup, and 
spurious polarization signal, can be reduced or avoided in interferometry. Being ground- 
based, it is controllable, and it can track the sky for an extensive period of time, as practiced 
successfully by the DASI team in measuring the CMB E polarization |[5j. Observational 
strategies of CMB interferometry experiments such as DASI, CBI, VSA, and AMiBA can 
be found in Ref. jl^l and references therein. 

In contrast to single-dish experiments which measure or differentiate the signals in indi- 
vidual dishes, an interferometer measures the correlation of the signals from different pairs 
of the array elements. The correlation output is called the complex visibility. In typical 
interferometric measurements, the field of view is small so that the sky can be treated as 
fiat. Therefore, the complex visibility is simply the two-dimensional Fourier transform of 
the intensity fluctuations on the sky convolved with the primary beam of the interferometer. 
The capability of directly sampling the Fourier modes allows for a simple estimation of the 
anisotropy power spectrum from visibility data j2l| and an efficient separation of the E and 
B polarization power spectra i^l- However, as the sky coverage is increased in mosaicking 



observations which combine several contiguous pointings of the telescope, the curvature of 
the sky becomes significant even though the sky is locally fiat at each pointing. Here we 
apply the results in the previous sections to extend the fiat-sky formalism of the complex 
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visibility to including the curvature effect of the sky. For further discussions about large- 
angular-scale interferometry in which three-dimensional Fourier transforms are involved, the 
reader may refer to Ref. 

Let us consider a two-element interferometer and a monochromatic electromagnetic 
source. The complex visibility is the time-averaged correlation of the electric field mea- 



sured by the two separated antennae pointing in the same direction to the sky j2J]. Usually, 
each antenna has a pair of feeds which are sensitive to orthogonal circular or linear po- 
larizations. For instance, if the dual-polarization feeds measure the right and left circular 
polarizations, then the output will be the four correlations (RR*), (RL*), (LR*), and (LL*). 
They can be related to Stokes parameters (T, Q, U, V). Denoting their associated visibility 
functions by (Vr, Vq, Vu, Vy), we have 

{RR*) = Vt + Vv, 

{LL*) = Vt- VV, 

{RL*) =Vq + iVu = 

{LR*) = Vq-iVu = V+2, (31) 

where we have neglected the parallactic angle of the feed with respect to the sky and the 
leakage from one polarization channel to the other polarization channel. The visibility 
functions are given by 

K(e, u) = ^j de'A{e'; e)X,(e')e2™-', (32) 

where u is the separation vector (baseline) of the two antennae measured in units of the 
observation wavelength, A denotes the primary beam with the phase tracking center pointing 
along the direction e, and Xg is the CMB sky. In Eq. (|32j) . dB^/dT is a conversion factor 
from the CMB temperature fluctuation to the brightness fluctuation given by 

^ 99.27- -rJy sr'i fxK-\ where x ~ 1.76 — , (33) 

dT (e^ - 1)2 ^ ' VlOOGHzy ' ^ ^ 

where B^, is the Planck function of the photon frequency u. 

Therefore, the response function of the interferometer is a complex function 

R\e'; e, u) = A{e'; 6)6^™-'^'. (34) 

Usually, the observation wavelength is much smaller than the size of the primary beam, 
dictating a small field of view. As such, for a single pointing, we can make the fiat-sky 
approximation by decomposing 

e' = e + /3, with ■ e = 0, and /3 < 1. (35) 

Here we assume a Gaussian primary beam given by 

A{(3,a)=exp(^-^y (36) 

Hence, by writing u ■ e' = u ■ (3 = u(3 cos a in Eq. (jSH), the interferometric response function 
becomes 

R^^{[], a) = A{f3, a)e*2™^=°^°. (37) 
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Note that in Eq. ()36|) we have adopted the interferometry convention that the primary beam 
does not carry the normahzation factor 1/(2710"^). Moreover, we have assumed a symmetric 
primary beam because a small beam asymmetry introduces only higher-order corrections 
which can be neglected as long as the length of the baseline is much bigger than the size of 
the dish. 

Similar to the previous case, substituting the beam (jHTj) in Eq. (fT^. making the approx- 
imation (|24j) in the integral, and using Eqs. (3.915.2) and (6.633.2) of Ref. jl8[ to perform 
the a-integration and /5-integration respectively, we obtain 



■s+m 



s+m 



(38) 



where lu = ^iru is the peak location of the primary beam. When the baseline length is much 
bigger than the antenna size, we have the limiting form 



L+m{llu(^b) - e"""VV27r//„(T2 for U^al > 1. 
Therefore, for long baselines, Eq. can be approximated as 



s+m 
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exp 



-a2(/-yV2 



which is actually a Gaussian beam centered at / = = 27cu. From Eq 
metric Gaussian window function is given by 

I 

tIG ^_2_4 r _2/,2 , ;2 ' ■ 



(39) 
(40) 

the interfero- 



47rV^xp -atil' + C 



^ 47rV,Sxp -al{l^ + ll] , . 

n=l n=l—l 

where z = lluO^l- Again, making use of Eq. (j^ . we find that for z < 1, 



n=l 

1-2 1+2 
II{Z)+2Y.II{Z)+ E ^nW 



47rV,^exp -at{l' + li) /o(2/U,'). 



(41) 
(42) 

(43) 



VI. FLAT-SKY APPROXIMATION 

We have discussed the full-sky CMB temperature and polarization correlation functions 
for asymmetric small-scale beams. In the limit of small sky coverage, their limiting forms 
must coincide with the results that are obtained in the fiat-sky approximation. In this limit, 
the sky can be treated as fiat, being spanned by a two-dimensional vector r. Hence, the 
CMB anisotropy and polarization fields are given by the Fourier transforms 

T(r) = j rfur(u)e-2™", 
Q(r) ± iU{r) = Jdu [e{u) t iB{u)] e^'^'^^e-^", (44) 

where 0u is the phase in the Fourier space given by the direction angle of u and 

(r*(u)r(w)) = Syy{u)6{u-w) where Y = T,E,B, 

(f*(u)E(w)) = Ste{u)6{u-w). (45) 

The power spectrum S{u) defined in the u-plane can be related to the angular power spec- 
trum Ci defined on the sphere by l{l + 1)Ci/2ti ~ 27iu'^S{u), with I ~ 27iu. 
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A. Single-dish Elliptical Gaussian Beam 



With an elliptical Gaussian beam with orientation angle ip given by 



1 



27r (TajCTy 



exp 



2al 



2(j2 

y 



where r = (r cos^,r sin^), the single-dish measurement gives 



(46) 



(47) 



where Xq = T, X±2 = Q T iU, and r is a pointing position on the sky. Using the expan- 
sion ()44|) and the Fourier transform of the response function 

(48) 



R^G{r,^) = J dui?^^(u,^)e-2™, 
the measured anisotropics can be written as 

X^^P{r,4j) = I (iw^^^(-w,^/')X,(w)e^"'^-e-2-™-r, 



(49) 



where Xq = T and Xj-2 = E ±iB. Hence, using Eq. fl45|) the two-point correlation function 
with different orientation angles at each point is given by 

{XT''*{r',^')Xr''{r,^P)) 
= J dwwSs'siw) J 0we-'(^'-^)'^-e-2™l'-'-'-l™^^-i?^^*(-w,^')^''^(-w,^), (50) 

where Sqq = Stt, S±2±2 = See + Sbb, S±2^2 = See - Sbb, Sn+2 = Ste, and 0w is the 
angle between r — r' and w. Expanding (see Eq. (8.511.4) of Ref. [l8|) 



-2niw\r —r \ cos f 



J2 (-^)™J„.(2™|r'-r|)e- 



(51) 



m=— oo 



and using Eq. (3.915.2) of Ref. [18| to perform the (/(-integration, we finally obtain 
(Xr'^*(r',^')^r"(r,^)) 



dw wSs'siw) exp 



TT 



-1^ 



^_i(s'_^+m)(V''+V)/2 ^ 

(52) 



r^Jm {27TW |r' - r|) I(s'-s+m)/2 (27r^u;^(o-^ - al) cositp' - tp) 
At zero lag, the non-vanishing contribution is from m = since Jm(0) = Sm,o and thus 
{\Xrnr,n') = J dwwSM exp [-2n'w\al + a^] 27r/o {27cV{cxl - a^)) , (53) 

which reproduces the result as found in Eq. (jSHI), with the window function given by Eq. ()3()|l 
and the replacement / = 2TTW. By taking the limit ax — >■ dy and using I(s'-s+m)/2{0) = 
6-m,s'-s, we obtain the standard two-point correlation functions with a Gaussian window 

(X^P *(r')XrP(r)) = J dwwSs'siw) exp [-ATrVa^) 27r(-z)^'-V,,_, (27ru;|r' - r|) . (54) 
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B. Interferometric Gaussian Beam 



In the flat-sky approximation (jH^j) . the complex visibility ()32|) is reduced to the two- 
dimensional Fourier transform of the Stokes parameter multiplied by the primary beam, 

OB, 



Kfr, u) 



dT 



dr'Aiv' - r)X,(r')e 



.l\ „2niu-r 



(55) 



where u is the two-dimensional projection vector of the baseline between two dishes in the 
r-plane and r is a pointing position on the sky. Therefore, we have 



K(r,u) = -^e~ 



dwA{u - w)Xsiw)e'''^'"e- 



(56) 



Hence, the two-point correlation function with different baselines at each point is given by 



(y;(r',u')K(r,u)) 



dT 



(W) X 



Assuming a Gaussian primary beam (jHUjl and using Eqs. (3.937.1) and (3.937.2) of Ref. [l 
to perform the (^-integration, we obtain for s' > s 

2 



(K*(r',u')K(r,u)) 



dB, 
dT 



g~2.i(u'.r'-u.r) / dwwSs's{w) CXp \-2lT^al{u'^ + + 2w^) 



3„4 



57r cr, 



where 



(p + iqf + 
p2 + (a + qY 

a = 27rw|r' — r| 



{s'-s)/2 



(58) 



p = 47r^Cr^w(u'cOS</)u' + MCOS(/)u), 

q = — 47r^cr^w(-u' sin 0u' + M sin 0u)- (59) 

As |r' — r| — s> and we set r = without loss of generality, we obtain the single-pointing 
visibility covariance matrices 

2 



(y;(uOK(u)) 



OB, 
dT 



(60) 



where U = \U\e^'^^ = u'e^'^^' + Me*"^", explicitly given by 

1 

\U\ = m'^ + + 2u'u cos(0u/ - 0u) % 

(pu = tan^^ [(u sin 0u' + u sin 0u) / {u cos 0u' + u cos 0u)] • 



(61) 



In a special case when |u'| = |u| = u, we have \U\ = 2m|cos[(0u' — 0u)/2]| and (pu 
(0u' + 0u)/2. Finally, the root-mean-square visibility is given by 

2 



VJr, u) 



OB, 
dT 



h{M'atuw), (62) 



which reproduces the result as found in Eq. (j25j) . with the window function given by Eq. 
and the replacement / = 2tiw. 
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VII. CONCLUSIONS 



Next-generation CMB experiments with high resolution and sensitivity will need to take 
into account the effect of beam asymmetry in interpreting observational data, while the 
pipeline for data analysis in CMB interferometric measurements using the fiat-sky approx- 
imation will have to be modified to including the effect of the sky curvature to deal with 
future observational data with large sky coverage. We have thus studied full-sky CMB 
anisotropy and polarization measurements with asymmetric beams. The measured Stokes 
parameters are well defined globally on the group manifold of the three-dimensional rotation 
and we have derived their correlation functions on the group manifold. The correlation func- 
tions are useful for constructing the full covariance matrices in the maximum-likelihood data 
analysis in large-sky CMB experiments with asymmetric window functions, especially when 
the beam is highly asymmetric and it is difficult to do beam symmetrization. Moreover, 
the domain of the CMB observable field is extended from the celestial sphere on which the 
field is expanded in terms of spherical harmonics to the rotation group manifold on which 
the field is expanded in terms of Wigner D-functions. As such, unbiased CMB angular 
anisotropy and polarization power spectra can be directly deconvoluted from observational 
data by converting the Wigner D-functions on the rotation group manifold. Work along this 
line is in progress. 
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